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Outline 

LET Flux Comparison - Interplanetary vs LEO at 28.5 deg & 51.6 deg 
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The Spacecraft Microelectronics Reliability and 
Qualification (MRQ) Process 

♦ Role of the mission and spacecraft specific single event effects 
(SEE) environment 

♦ Role of nuclear reaction and transport processes in 
determining the spacecraft SEE environment 

Natural vs. induced SEE environments ± nuclear reaction in 
shielding mass and avionics components 

♦ Why FLUKA (FLUktuierende KAskade)? 

Comparing in-flight single event upset rates with rates 
calculated using FLUKA LET spectra 

♦ Low-Earth Orbit (LEO) 

Space Shuttle and International Space Station (ISS) 
Calculated vs. observed single event rates 
^ LET spectra as a function of shielding mass 

♦ Geosynchronous and Interplanetary orbits 

Cassini, SOHO, Mercury Messenger, ETS-V, Thyraya 
^ Calculated vs. observed single event rates 
LET spectra as a function of shielding mass 

Calculating SEE environments and rates with FLUKA 

♦ Calculating LET Spectra 

♦ Calculating SEE rates 

♦ LET spectrum and SEE rate characteristics 

^ Shielding mass effects 

Contribution of various CR elements to the expected SEU 
rates 

The effects of high Z elements in microelectronic devices 

Discussion, Summary and Conclusions 



The Spacecraft Microelectronics 
Reliability and Qualification 
( MRQ ) Process 

Step 1- Mission Specific SEE Environment Specification 

Step 2 - Microelectronic device SEE Characterization 

Step 3 - Calculate (estimate) expected in-flight device SEE rates 

Step 4 ± Is the expected SEE rate acceptable in light of system 
safety and reliability requirements? 

How do we know all this works (method validation)? 

Why the interest in FLUKA? 


The Spacecraft Microelectronics Reliability and 
Qualification (MRQ) Process 


Step 1- Mission Specific SEE Environment Specification 

♦ Natural SEE Environment 

Galactic Cosmic Ray (GCR) and Trapped Radiation environment definitions from the 
U.S. Naval Research Laboratory CREME 96 web page, (https://creme96.nrl.navy.mil/) 
(for solar minimum) 

♦ Spacecraft Specific Induced SEE Environment 

Results when the natural SEE environment is processed by: 

4 1 ) spacecraft structural, consumable, and shielding materials 
4 2) avionics device materials 

tsr Calculated using a nuclear reaction and transport model ± FLUKA in this case 

The result is a differential, f[LET], or integral F[LET] probability distribution function 
providing the particle flux as a function of particle LET - isotropic to first approximation 

LET units are in (MeV cm 2 )/mg (Si) 

Step 2 - Microelectronic device SEE Characterization 

♦ Device (chip) level heavy ion and/or proton testing 

^ Cross section for single events vs. charged particle Linear Energy Transfer (LET) in the target device 
at one of several accelerator facilities 


The Spacecraft Microelectronics Reliability and 
Qualification (MRQ) Process 

Step 3 - Calculate (estimate) expected (in-flight) device single event rates 
using: 

♦ 1 ) Differential form of the SEE environment definition, f[LET], solar minimum 

f[LET] is isotropic (ICRU definition) 

Scale the simulation particle flux/fluence to on-orbit flux/fluence 

♦ 2) The integral form of the device directional cross section 


WHY FLUKA? 


• Microelectronics are evolving beyond the limits of established single event 
test and verification methods 

♦ The size of microcircuit elements continue to shrink 

♦ Microcircuit element packing density and VLSI circuit complexity continue to increase 

♦ High atomic number elements like Ni, Cu, W and Hf are introduced over a range of size 
scales introducing nuclear reaction issues 

♦ Many of the assumptions underlying affordable versions of the established test methods 
(JEDEC EIA/JESD57 and/or ASTM F1 192 -00(2006)) and associated SEE rate models 
are increasingly unreliable or not applicable 

♦ (Johnston, 1998), (Shaneyfelt, Schwank, Dodd, Felix, 2008), (Lacoe, 2008), (Schrimpf, Warren, Ball, Weller, Reed, Fleetwood, 

Massengill, Mendenhall, Rashkeev, Pantelides, Alles, 2008), (Muntenau, Autran, 2008), (Fulkerson, Nelson, Carlson, 2006), ( 

Reed, Kinnison, Pickel, Buchneer, Marshall, Kniffin, LaBel, 2003) 

• Increasing community interest in Monte Carlo nuclear reaction and transport 
codes to understand and evaluate SEE processes in modern microelectronic 
devices 

♦ Spacecraft and planetary/lunar/asteroid shielding mass effects 

♦ Energy/charge deposition 

♦ Nuclear reactions in or near microcircuit elements 

♦ Multi-node charge collection 

♦ Simulation of heavy ion testing 

♦ Multiple bit upset analysis 

♦ Simulation of low energy and high energy proton testing 

♦ □)RVWHUDD2i|iHLOOan.RXEDnnnnnnnnnn6FKULPSlnn:DUUHQnn%DOOnn:HOOHUnn5HHGnn)OHHWZRRGnnODVVHQ 
Pantelides, Alles, 2008), (Warren, Sierawski, Reed, Weller, Carmichael, Lesea, Mendenhall, Dodd, Schrimpf, Massengill, 

Hoang, Wan, De Jong, Padovani , Fabula, 2007) (Skutnik, Laioie, 2006) 



Why FLUKA? - Secondary Particle Production in Spacecraft 
Shielding , Avionics Components , and the Atmosphere 


• Three basic processes 

♦ Energy loss (dE/dx) by ionization of material along the particle 

track 

«■ Direct ionization effects ± linear energy transfer (LET) 

♦ High energy collision (inelastic) with spacecraft materials nuclei 

Nuclear reactions initiate secondary particle showers 
Proton and neutron SEE effects are often the result of direct 

nuclear reactions 


4 reactions in or near the device sensitive volume 
Further collisions of secondary particles with spacecraft nuclei 
leading to expansion and propagation of the secondary particle 

shower 

Secondary particles can produce direct ionization and more 

nuclear reactions 

4 Recoil nuclei have short range but high LET 

♦ Collisions with material nuclei produce displacement damage 

• Basic physics is similar to cosmic ray air showers except ± 

♦ Higher atomic number spacecraft target nuclei produce more 
secondary particles 

♦ Density of spacecraft materials » density of air 

• Minimum and Median Shielding Mass 

♦ Inside ISS - US Lab 
Minimum = 10 g/cm 2 Al 
Median = 40 g/cm 2 Al 

♦ (DUWK^VHDWPRVSKHUH DW □ ni NP DOWLWXGH □ 
Minimum = 56 g/cm 2 air 
Median = 84 g/cm 2 air 


This centra i Au on Au collision at 4 GeV/nucleon 
kinetic energy comes from the January 1996 run of 
E895. ionization clusters in the EOS Time Projection 
Chamber are shown in red. The chamber iength is 
about 2 meters. 

Analysis of Data from Experiment E895 at the Alternating 
Gradient Synchrotron; Dieter Best, Gulshan Rai, and Hans - 
Georg Ritter, Lawrence Berkeley National Laboratory, 1998 




WHY FLUKA? - Demonstrated successful applications include: 


♦ Cosmic ray physics 

♦ Neutrino physics 

♦ $FFHOHUDWRUO GHVLJQU □ 


Spacecraft SEU rate calculations and 
comparison with in-flight rates 


Success Metric Ratio = In-Flight SEU Rate/Predicted SEU Rate 

Success Criteria = 0.1 < (In-Flight Rate/Predicted Rate) <10 

This section presents results only ± details of calculation methods, 
approximations used, and results will be treated in the section 3 
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Weekly Average SEU Count by MDM Group 
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Enhanced MDM 
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ISS & Shuttle Success Metric Plot for 2 different spacecraft ; 4 
different parts, 3 different shielding mass environments 



See Appendix 2: Tables 1-5 for case 
specific details 


FOM, figure of merit calculation 


FLUKA sphere target calculation 


FLUKA slab target calculation 


1) FLUKA and FOM meet 
the success criteria or 
are very close in all 

cases 

2) FLUKA results for the 
median shielding mass in 
the spherical target and 
the simple shielding 
mass in the slab target 
are in agreement 


Shielding mass effects: ISS FLUKA spherical 

target calculation 


Y axis - ISS MDM DRAM, 
upsets per week per MDM 


Tl (4Mx 4) 100 

Tl SMK416400 


FLUKA calculation 
methods predict shielding 
mass effects on SEU rate 
accurately for the standard 
ISS MDM DRAM but not as 
accurately for the 
Enhanced MDM DRAM. 


80 

i- 


Tl (1M x 4) \ 

TMS44400 \ 

60 - 

STDmdmratelOj 

>ee< 


Xaxis - Median ENHmdmrate20 i 

shielding mass in 
grams/cm 2 


Shielding mass 

Rate Ratio -(10 g/cm 2 

Rate)/ (40 g/cm 2 Rate) 

Device 

Rate Ratio - Flight 

Rate Ratio ± FLUKA 
Calculation 

Tl (1M x 4) 
TMS44400 

1.2 

1.2 

Tl (4M x 4) 

Tl SMK416400 

0.9 

1.8 





CEO/Interplanetary Success Metric Plot for 5 different 
spacecraft , 6 different parts and shielding mass environments 


Mercury Messenger 

Right Circular Cylinder (RCC) Target, T/W=1) 



RatioFOM is Rsphi 


Mercury messenger shielding mass and tungsten over layer effects: 


1x10 


.-10 


Y axis - Mercury Messenger 
SRAM upsets per bit day, 
RCC, T/W=1 


MMRCCrate j 

MMRCCrateWj 

-H-+ 


For Mercury Messenger SRAM, a 
FLUKA rate calculation, using an 
approximate sensitive volume 
geometry, predicts a very different 
SEU rate dependence on shielding 
mass when the 1 micron Tungsten 
over-layer is present above the Si 
detector shell 
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X axis ± g/cm 2 Al median shielding mass 
(spherical target , cosine and solid angle corrections) 


Calculating SEE Environments and 

rates with FLUKA: 


Calculating LET Spectra 

Calculating SEE rates 

LET spectrum and SEE rate 
characteristics 


FLUKA- Comprehensive basic physics with experimental 

validation 
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• FLUKA is a multipurpose Monte Carlo energetic particle interaction and transport 

code 

♦ Monte Carlo code (more than 500,000 lines of Fortran Code) 

♦ Theory driven/experimentally benchmarked ± well tested physics based microscopic 
models (not semi-empirical look-up tables) 

Nucleus-nucleus collisions and secondary particle production included explicitly 

♦ Benchmarked/verified extensively with high energy accelerator data 

not yet benchmarked for spacecraft SEE/TID/DD processes of interest 
FLUKA 2006.3b only partly successful (MRQW 2007) 

♦ FLUKA is not a tool kit ± the physical models are fully integrated 

♦ Full development of hadronic (secondary particle) showers 

♦ Complex user generated target geometries and target materials are possible 

• Basic References 

♦ FLUKA version ± Fluka2008.3b results reported here 
FLUKA results reported at MRQW 2007 generated with an earlier version of FLUKA 

♦ □ 7KHD)/8.$\3FRGHD □ 'HVFULSWLRQr\DQGaEHQFKPDUNLQJ 3 D*nD %D WWL VWRQL □ □ 6 □ □ OXUDUR □ C 
Ferrari, S. Roesler, A. Fasso', J. Ranft; Proceedings of the Hadronic Shower Simulation Workshop 2006, 
Fermilab 6-8 September 2006, M. Albrow, R. Raja eds., 

AIP Conference Proceeding 896, 31-49, (2007) 

♦ "FLUKA: a multi-SDUWLFOHDWUDQVSRUWDFRGH 3 D$nn)DVVRCnD$nn)HUUDULnn-nD5DQIWnnDQ 
CERN-2005-1 0(2005), INFN/TC_05/11, SLAC-R-773 

♦ http://www.fluka.org/fluka.php 
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“System Requirements ” 

• Run Times on the JSC/ES4 blade server 

♦ For a given CR spectrum and shielding mass or target configuration: 

We run each cosmic ray element as a single job 
•»“ Using a blade server we can submit all 24 to 21 jobs at once and they will execute in 

parallel 

4 :HUKDYHU □ □ GXDOU TXDGUFRUHU 3 EODGHV'U GHGLFDWHGD WRU )/8 . $□ UXQVl 
4 (DFKn 3 FRUH'DLVnHTXLYDOHQWn WR UDQnLQGHSHQGHQWn SURFHVVRUn IR 
4 So, we have 24 processors available and can run up to 24 jobs in parallel 

4 All the runs ( 24 to 27 batch jobs submitted in parallel) generally complete wjtjTin 
one week gf_submissign (see chart 18) 

Fora fixed number of primary particles, run time increases dramatically with particle 
atomic number 

• Without the parallel processing capabilities of a blade server, execution 
times are too long to be of any practical use in most cases 

♦ Unless only one CR or Trapped particle species is of interest 

• JSC/ES4 blade server specification: 

♦ Blade Hardware - HP ProLiant BL460c G1 x64 
2 Quad-Core Intel Xeon 2.6 GHz processors with 32 GB RAM and 146 GB storage 

♦ OS - RedHat ES5 version 5.2 
Batch job queue management - MCS Portable Batch System 
Fortran G77 compiler 


Calculating Spacecraft SEE Environments with FLUKA 


8VHnnJHQHULFDVODEnRUnFRQFHQWULFnVKHOOn 3 VSDFHFUDIW'nJHRPHWI 
VLOLFRQnnRUnRWKHUno 3 VFRULQJ'nGHWHFWRUVnEHWZHHQnOD\HUVDRIn 
shielding mass 

♦ Simple, well defined shielding mass distribution and median shielding mass for each 
detector shell in each geometry 

♦ For the slab target geometry FLUKA simply fires particles into the center of the slab at 
normal incidence 

♦ For the concentric shell target, FLUKA fires randomly directed particles, selected from the 
spacecraft natural SEE environment, into the spherical shell structure from the outside (an 
ICRU isotropic flux) 

FLUKA utilities randomly samples natural GCR and Trapped particle spectra 
to specify particle, particle kinetic energy, and particle direction 

♦ For each particle, FLUKA calculates through the target structure along the particle 
track: 

Energy loss (LET) of primary particles 

Nuclear reactions and reaction products (secondary particle showers) 

Energy loss (LET) and further nuclear reactions of secondary particles 

A FLUKA utility generates the final product, i.e. the LET spectrum entering 
each Si detector shell which includes all contributions to LET from both 
primary particles and secondary particles formed in shielding 

♦ FLUKA also calculates nuclear reactions and recoil products interior to the Si layer 
but efficient methods for scoring those contributions to the detector shell LET 
spectrum are still in development at this time 

♦ FLUKA reports both forward and backward going particles with respect to the 
primary beam direction ± only forward going are reproted here. 


Calculating Spacecraft SEE Environments with FLUKA 

As employed here, FLUKA calculates an estimate of the LET spectrum 
internal to a microelectronic device die material located internal to 
spacecraft shielding mass 

♦ Device or system SEE rates are then calculated using: 

«■ the FLUKA detector shell LET spectrum 
device SEE characterization data 
directional cross section models 

♦ We do not calculate SEE rates during the FLUKA calculation based on specific sensitive 
volumes imbedded in the Si detector shells at this time 

The usual assumption is made 

♦ Energy deposition and charge production in the target is proportional to the product of 
ion LET and ion track length and LET is assumed constant 

:KDWflVnQHZnDQGDGLIIHUHQWnKHUH ,, D 

♦ Detailed treatment of nucleus-nucleus collisions in spacecraft SEE transport calculations 
with comprehensive treatment of all secondary particles 

♦ As employed here, FLUKA calculates an estimate of the LET spectrum internal to a 
microelectronic device die material located internal to a spacecraft 

Nuclear reactions internal to the detector shell are counted but the contributions of 
the nuclear reaction products to the detector shell LET spectra are not at this time 

We compare the FLUKA SEE rate with the in-flight SEE rate for the 
applicable spacecraft median shielding mass 

♦ Success Criteria - As a minimum, the on-orbit SEE rate calculation method should 
provide SEE rate estimates accurate to within a factor of 1 0 at one standard deviation 
when compared to available in-flight data 


0.1 < (In-Flight Rate/Predicted Rate) <10 
Success Metric Ratio = In-Flight Rate/Predicted Rate 


FLUKA At Slab Target 

10 micron silicon detector layers 


Particle Beam 

* 



X axis - distance from center of cylindrical AL slab target in cm 
(The actual slab target diameter used was 100 cm) 

Y axis - distance form top surface of slab to Si detector slab in cm 


FLUKA Concentric Ai Sphere Target in Cross Section 

100 micron Si detector shells, polar coordinates 


o 

5000.0 
5018.9 

5030.0 
5033.7 
5035.6 

5037.1 
50373 
5037.4 



The volume of the sphere at radii smaller then 
5000 cm is treated as a perfect particle absorber 
in all FLUKA runs 


FLUKA 2008.3b Calculation Details 


Number and identity of GCR primary particles 

♦ P, He, C, O, Mg, Si, Fe, and Zn only, no higher Z elements 

♦ Number of primaries particles for FLUKA runs are always greater than the expected 
CREME-96 weekly fluence (#/cm 2 week) to assure adequate and physically realistic 
statistics 


Element 

H trapped 

H 

He 

C 

0 

Mg 

Si 

Fe 

Zn 

CREME fluence 
ISS orbit 

2.5 xIO 7 

3.5 x10 s 

4.7 x 10* 

1.3 x10 s 

1.3 x10 s 

261 

193 

157 

0.14 

CREME fluence 
Interplanetary 

0.0 

2.82 x10 s 

2.7 x 10 s 

6.9 x10 s 

7.1x10 s 

1.4 x10 s 

996 

731 

0.64 

FLUKA primaries 
ISS orbit 

9x 10 7 

9 x10 s 

9 x10 s 

9 x10 s 

6.8 x10 s 

4.5 x10 s 

4.5 x10 s 

2.3 x10 s 

9x10* 

FLUKA primaries 
Interplanetary 

0.0 

9 x10 s 

9 x10 s 

9 x10 s 

6.8 x10 s 

4.5 x10 s 

4.5 x10 s 

2.3 x10 s 

9x10* 


Slab Target - particles flux at normal incidence - slab 10 meters in diameter 

Spherical Shell Target ± FLUKA isotropic (ICRU) flux utility provides an isotropic flux to the 
exterior surface of the sphere 

♦ Inner radius of concentric spherical shell structure = 5,000 cm = 50 meters 

LET scoring ± 0.0010 or 0.0100 cm thick Si scoring targets for both slab and concentric 
spherical shell targets for each shielding mass thickness 

♦ LET from all particles and interactions to include recoil products from proton and neutron 
reactions in the Si scoring target 


FLUKA 2008.3b Calculation Details 


• FLUKA Target Shielding mass for each Si detector layer 

♦ From outside to inside in g/cm 2 total shielding mass 
(overlying Al or PE + Si scoring target mass) exterior to 
each Si scoring detector (SiDet) 

♦ Spherical shell and slab targets 

♦ Total areal shielding mass - from outside to inside with 
respect to the entering particle beam 


FLUKA Target 

SiDetl 

SiDet2 

SiDet3 

SiDet4 

SiDet5 

SiDet6 

SiDet7 

SiDet8 

Slab 

0.1 

0.5 

1.0 

5.0 

10.0 

20.0 

50.0 

100 

Sphere (minimum) 

0.1 

0.5 

1.0 

5.0 

10.0 

20.0 

50.0 

100 

Sphere - median (cosine 

0.14 

0.70 

1.40 

6.90 

13.7 

27.3 

68.1 

137.2 

correction only) 









Sphere ± median 

0.15 

0.81 

1.6 

7.9 

15.6 

31.1 

77.5 

156.2 


cosine and solid angle 
corrections 











FLUKA differential LET spectra - ISS Orbit 
SLAB Target 


Y axis: on-orbit 
particle flux to each Si 
detector in units of#/ 
[cm 2 /week] per LET 
unit ± Shielding 
values per page 27 for 
SLAB target 


SiDetl 

SiDet2 

SiDet3 

SiDet4 

SiDet5 

SiDet6 

SiDet7 

SiDet8 



(LETmeanBIN) 


X axis: particle LET in MeV x mg/(cm 2 ) (Si) 


FLUKA differential LET spectra - ISS Orbit - SPHERE 
Target 


Y axis: on-orbit particle flux 
to each Si detector in units 
of# / [cm 2 /week] per LET 
unit ± Shielding values per 
page 27 for SPHERE target; 
median mass, cosine and 
solid angle corrections 

SiDetl 

SiDet2 

SiDet3 

SiDet4 

SiDetS 

SiDet6 

SiDet7 

SiDet8 



(LETmean BIN) 


X axis: particle LET in MeV x mg/(cm 2 ) 


28 

Device Directional Cross Sections 


• Problem - isotropic flux on a generally anisotropic target 


Generalized Rate Equation 

♦ Polar angle ( 


ISS orbit SEU rates for a simple isotropic Step Function o 

SLAB SPHERE 




X axis ± g/cm 2 Ai median shielding mass 

(For the spherical target , cosine and solid angle corrections per page 27) 
Y axis - number of step function isotropic target hits per week; 


Which energetic nuclei contribute most to the step 
function SEU rates in LEO ? 


Alsphrate^ 

Alsphratet) 


Alsphratei} 

“ 0 " 

Alsphrate4 



AknedianCO ®vf \ 


AltredianCOSA'f, 


AlmedianCO ©vfj 


X axis ± g/cm 2 Al shielding mass, spherical target with cosine and view factor corrections per page 27 
Y axis - number of step function isotropic target hits per week; 


FLUKA differential LET spectra - GEO/ Interplanetary - 
SPHERE Target - Silicon detector shells only 


Y axis: on-orbit particle flux 
to each Si detector in units 
of#/ [cm 2 /week] per LET 
unit ± Shielding values per 
page 27 for SPHERE target; 
median mass, cosine and 
solid angle corrections 

SiDetl 

SiDet2 

SiDet3 

SiDet4 

SiDetS 

SiDet6 

SiDet7 

SiDet8 



X axis: particle LET in MeV x mg/(cm 2 ) (Si) 


FLUKA differential LET spectra - GEO/inter planet ary - SPHERE 
Target - Silicon detector shells with i micron tungsten over - 
layer 


Y axis: on-orbit particle flux 
to each Si detector in units 
of# / [cm 2 /week] per LET 
unit ± Shielding values per 
page 27 for SPHERE target; 
median mass, cosine and 
solid angle corrections 

SiDetl 
SiDet2 
SiDet3 
SiDet4 
SiDet5 
SiDet6 
SiDet7 
SiDet8 



(LETmeanBIN) 


X axis: particle LET in MeV x mg/(cm 2 ) 



CE O/Interplanetary SEU rates for a simple isotropic Step 
Function o 
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X axis ± g/cm 2 median Al shielding mass , spherical target with cosine and solid angle corrections per page 27 
Y axis - number of isotropic target hits per week; 


Which Cosmic ray nuclei contribute most to the step 
function SEU rates ? 


Full FLUKA Free Space GCR set 



Removing Free Space GCR Fe 



Removing Free Space GCR Fe, He, P 



Alnc(hanCOS>vf i 


X axis ± g/cm 2 median Al shielding mass, spherical target with cosine and solid angle corrections per page 27 
Y axis - number of isotropic target hits per week; 


The answer depends on the step function threshold and 
the shielding mass. 


Discussion Summary and Conclusions 


The FLUKA LET spectra in combination with heavy ion test data 
successfully predicts on-orbit SEE rates for several different CMOS 
SRAM, DRAM, and SDRAM components 

♦ Success criteria 1 ± The FLUKA predicted rate is within a factor of 10 of 
the in-flight rate in many cases using approximate treatments of the 
directional cross section 

Success criteria 2 - Shielding mass effectiveness may still be overestimated in some 
cases; accurate for the standard ISS MDM CMOS DRAM 

♦ Compares favorably with the well documented Figure of Merit method for 
those cases where the Figure of Merit is applicable 

♦ Predicts an increase in SEU rates produced by inclusion of a high Z 
element (W) in micro-device structure , all else being equal 

♦ Predicts important effects of the shielding mass on the SEU rate increase 
produced by a high Z element (W) in the micro-device architecture 

So what happened between MRQW 2007 and MRQW 2009? 

♦ FLUKA 2006.3b vs. FLUKA 2008.3b 

♦ 7KHU )/8.$UFRQVRUWLXPnXSGDWHGU WKHD 3 86(pt$XWLOLW\n WRUSURYL 
more accurate/complete reporting of high LET particles (personal 
communication, Professor L. S. Pinsky, Physics Department, University of 
Houston, FLUKA consortium) 


Discussion Summary and Conclusions 

• The FLUKA nuclear transport and reaction code can be developed into a 
practical tool for calculation of spacecraft and planetary surface asset SEE 
and TID environments 

♦ Nuclear reactions and secondary particle shower effects can be estimated with 
acceptable accuracy both in-flight and in test 

More detailed electronic device and/or spacecraft geometries than are reported here 
are possible using standard FLUKA geometry utilities 

Spacecraft structure and shielding mass 

*»" Effects of high Z elements in microelectronic structure as reported previously (Schrimpf, 

Warren, Ball, Weller, Reed, Fleetwood, Massengill, Mendenhall, Rashkeev, Pantelides, Alles, 2008), (Warren, 
Sierawski, Reed, Weller, Carmichael, Lesea, Mendenhall, Dodd, Schrimpf, Massengill, Hoang, Wan, De Jong, 
Padovani , Fabula, 2007) 

♦ Median shielding mass in a generic slab or concentric sphere target geometry are 
at least approximately applicable to more complex spacecraft shapes 

Need the spacecraft shielding mass distribution function applicable to the 
microelectronic system of interest 

♦ SEE environment effects can be calculated for a wide range of spacecraft and 
microelectronic materials with complete nuclear physics 

^ Evaluate benefits of low Z shielding mass can be evaluated relative to aluminum 
Evaluate effects of high Z elements as constituents of microelectronic devices 

• The principal limitation on the accuracy of the FLUKA based method 
reported here are found in the limited accuracy and incomplete character of 
affordable heavy ion test data 

♦ To support accurate rate estimates with any calculation method, the aspect ratio of 
the sensitive volume(s) and the 


Appendices 




Appendix 1: FLUKA based in-flight SEE 

rate calculation details 


FLUKA Physics Overview 

• Heavy ion interactions models 

♦ E> 5 GeV/n 

Dual Parton Model, DPMJET-III 

♦ 0.1 GeV/n < E< 5GeV/n 

Relativistic Quantum Molecular Dynamics, RQMD 2.4 

♦ E< 0.1 GeV/n 

i3r Boltzmann Master Equation (BME) Theory 
BME code ± not utilized in this study 

• Hadron - Nucleus Interactions 

♦ Resonance production and decay below a few GeV energy 

♦ Dual Parton model above a few GeV energy 

The PEANUT model includes a detailed Generalized Intra-Nuclear Cascade (GINC) 
and a pre-equilibrium stage 

Gribov- Glauber multiple collision model included in a less sophisticated GINC 

• Transport of charged particles in matter 

♦ Bethe-Bloch theory 

S/7e// and other low energy corrections, 

Density effects according to Sternheimer 

WKHRU\URlUPXOWLSOHU &RXORPED VFDWWHULQJ 
Restricted fluctuations (Landau fluctuations) 

Delta ray production and transport optional 


Typical FLUKA input (run) Physics Card Choices 


With the PRECISION default we request the following 
nuclear physics model 

1) Dual Parton Model Jet (DPMJET) 

with Relativistic Quantum Molecular Dynamics (RQMD) 

2) Evaporation , Coalescence , and Electromagnetic Dissociation 
are enabled 

3) The Peanut model is activated at all energies 


TITLE 






CREME 96 GCR Spectrum 




DEFAULTS 





PRECIS 10 

BEAM 

-100.0 




HEAVYION 

BEAMPOS 

0.0 

0.0 5037.2 




HI-PROPE 

26.0 

56.0 0.0 0.0 

0.0 

0.0 

EVENTYPE 


2.0 



DPMJET 

DPMJET 

0.0 

0.0 0.0 



my card 

PHYSICS 

3.0 

0.0 0.0 0.0 

0.0 

0 . 

EVAPORAT 

PHYSICS 

1.0 




COALESCE 

PHYSICS 

2.0 

0.0 0.0 0.0 

0.0 

0 . 

EM-DISSO 

PHYSICS 

1000. 

1000. 1000. 1000. 

1000. 

1000. PEATHRES 


SEE Rate Calculations 


Isotropic Target, (Edmonds, Barnes and Scheick, 2000) 

♦ 


SEE Rate Calculations 

Right Circular Cylinder (RCC) Target 

♦ Barak, Akkerman, 2005 

♦ Akkerman, Barak, 2002 

♦ Barak, 2001 

♦ Barak, Reed, LaBell, 1999 

Note that we use the average (first moment) cord length for a given 


Appendix 2: Spacecraft, Device Data, 
ED AC Protocols, and Shielding Mass 

Distribution Functions 


iSS DRAM and SDRAM Characteristics 
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□ ISS standard multiplexer de-multiplexer (std-MDM) input-output control 
unit (IOCU) board 

□ Texas Instruments (Tl) (1M x 4) 4 Mbit CMOS DRAM TMS44400 

□ 3.3554 x 10 7 bits in each std-MDM 

□ EPIC Process 

□ 0.9 micron device scale 

□ 5 V 

□ ISS enhanced MDM DRAM (enh-MDM) IOCU board 

□ Tl (4Mx4) CMOS DRAM Tl SMJ416400 

□ 1 .342 x 10 8 bits in each enh-MDM 

□ 5V 

□ ISS High Rate Communications Outage Recorder Samsung KM44S32030T- 
GL CMOS SDRAM 128 Mbit (Rev. A, 34M/4) 

□ ISS High Rate Communications Outage Recorder (HCOR) 

□ 2.1 15 x 10 11 bits SDRAM in one HCOR (constant since launch! losses made 
up from reserves) 

□ 0.35 micron process/device scale 

□ 3.3 V 
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fSS DRAM and SDRAM Weibull Parameters 


STD MDM DRAM cross section vs LET 

lostd:=0.99 astdsat :=3.0010 7 wstd:=7.7 zstd:=1.3 


astdMDM (1) : = astdsat ■ 

ENH (Enhanced) MDM rates - SM J416400 

loenh:=0.42 aenhsat : = 1.1010 




( 1-0.1- lostd 

1 - exp 

— 

: — 1 i 



^ wstd J _J 


8 


wenh : = 0 . 8 zenh := 1.7 


aenhMDM (1) := aenhsat 


1 - exp 


( 1-0.1- loenh \ 


wenh 


J 


enh "| 

111 


STD MDM DRAM (TMS44100DM-80) characterization - Harboe 
Sorensen , Muller, Daly, Nickson, Schmitt, Rombeck 1991) 


Heavy ion data for Samsung 128Mbit (Rev. A, 34M/4) S DRAW: average values 

from Henson, MacDonald, and Stapor, NSREC Workshop 1999 - Samsung KM44S32030T-G 

Fig. 6, static tests 


HCOR SDRAM upper bound cross section vs LET 1:= 1 .. 500 

— 8 

losdraml := 13 asatsdram! ;= 5.85910 wsdraml := 30 zsdraml := 1 


aHCORhigl^l) : = asatsdram 1 


1 - exp 


1-0.1- losdraml \ zsdram j 1 
wsdraml J | 


HCOR SDRAM lower bound cross section vs LET 

— 8 

losdram2:=14 asatsdram 2 := 1 .5631 0 wsdram2 := 30 zsdram2:=l 


aHCORlo^l) := asatsdram 2' 


1 - exp 


10.1- losdram2Y sdrarri | "I 
wsdram2 J j 


Heavy ion data for Samsung 128Mbit (Rev. A, 34M/4) SDRiA&d lvalues from SEAKER repoi 
SDRL number MD005, SDS Number SS-EE-008, transmittal number 00-HCOR-086, PDC 
Number FM27257, 


ENH MDM DRAM characterization - Brown, R., IBM Manassas Test 
Report 2/23/93, and Falguere, Duzellier, Ecoffet, Tsourilo, 2000 

HCOR rates -HCOR SDRAM cross section vs LET (SEAKR- low ke ions, normal incidenc 
only except for one 40.5 degree angle and one 30.2 degree angle) - Boeing Dr. Tiru 
Rao-Sahib, and Aerospace Corp. Dr. Rocky Koga; U.C. Berkeley Cyclotron facility August 
2 , 2000 . 


- Q 

losdram3:= 1.95 asatsdram3 := 1.86310 wsdram3 := 30 


aHC OR3(l) : = asatsdram 3 ■ 


1 - exp 


1-0.1- losdram3 Y 


sdram3 ] 


wsdram3 j j 


zsdram3 := 1.9 


ISS: Memory device information and Error 
Detection and Correction (ED AC) Protocols 


• Standard MDM Memory Information (MDM-4, MDM-10, MDM-16) 

♦ Part Number = TMS44400; Part Type = 1Mx4 DRAM 

♦ Manufacturer = Texas Instrument (Tl) 

♦ Device count (per MDM) = 8 

♦ Total Bits (per MDM) = 33,554,432 

♦ Memory Scrub Rate = memory initiated every 8 usee, 8.2 seconds total to scrub entire DRAM memory 

• Enhanced MDM Memory Information 

♦ Part Number = SMJ416400 ;Part Type = 4Mx4 DRAM 

♦ Manufacturer = Texas Instrument (Tl) 

♦ Device count (per MDM) = 8 

♦ Total Bits (per MDM) = 134,217,728 

♦ Memory Scrub Rate = initiated every tbd usee, 28.6 seconds total to scrub entire DRAM memory 

• HCOR Memory Device Information 

♦ Part Number = KM44S32030A T 

♦ Part Type = 128MB SDRAM 

♦ Manufacturer = Samsung Memory 

♦ Array = 140 128Mbit devices (70 stacks of 2 high) 

♦ Memory Utilization = 211,452 Mbits (this is the amount of scrubbed memory available for users, 
accounting for spare boards, image RAM, etc.) 

♦ Scrub Rate = 3.25 usec/row, 7.3 min/board, 96 minutes total to scrub entire memory (assumes minimal 
errors encountered. One scrub in progress at any given time). It is not possible to generate position data 
with the HCOR SDRAM error counter data because of the delay (up to 95 minutes) between an SEU event 
occurring and the ED AC scrub correcting and counting the error. Refresh rate for error counter telemetry 
is Ihz prior to January 12, 2005 and 0.1 hz from January 12, 2005 to the present. 


Fraction of solid angle< T 
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ISS Shielding Mass Distribution Functions: 

1) Multiplexer/De-multiplexer (MDM) intrenai vs external 

2) High Rate Communication Outage Recorder (HCOR) 



0.1 1 10 100 


US Lab Rack #3, HCOR location inner/outer boundry 



SMdmassI^JShldmassEXT 


Shielding mass in grams per square cm 

HCOR SDRAM Structural Shielding Distributions. 

Pendleton, G. N.; 3 5DGLDWLIBQse Sector Shielding Track Analysis Tool Results for the International Space 
6WDWLRQolsa Corporation Report Colsa-RTD-ISS-DR-03-008-DOC-B, Colsa Corporation, Huntsville Alabama, 
USA, March 2003 (Marshal Space Flight Center cooperative agreement NCC8-200, J. W. Watts Jr. Project 
Technical Representative) 


Equivalent Aluminum Thickness, T (g/cm2) 

MDM DRAM Structural Shielding Distributions. 


Table 1: Std-MDM DRAM and enh-MDM DRAM upset rates are in reasonable 
agreement with FOM predictions and FLUKA predictions 


Device 

Number of 

bits 

Per MDM 

Median 

Shielding 

Mass 

Device 

FOM 

FOM 

calculated 
upsets per 
week per 

MDM 

FLUKA Sphere Target 
RCC calculated upsets 
per week per MDM 

(SV aspect ratio)* 

Range of 
Observed In- 

Flight 

Upsets per week 

Tl (1M x 4) 
TMS44400 

33.55x10° 

40 g/cm 2 

1.93x10* 

16 

17 

(1) 

15 to 18 

Regression line 

Tl (1M x4) 
TMS44400 

33.55x10° 

10 g/cm 2 

1.93x10* 

59 

21 

(1) 

18 to 22 

Regression line 

Tl (4M x 4) Tl 
SMK416400 

1.342x10° 

40 g/cm 2 

8.98x10-° 

2 

26 

( 2 ) 

(Falguere, Duzellier, Ecoffet, 
Tsourilo, 2000), 

3 to 4 

Regression line 

Tl (4M x 4) Tl 
SMK416400 

1.342x10° 

10 g/cm 2 

8.98x10-° 

9 

48 

( 2 ) 

(Falguere, Duzellier, Ecoffet, 
Tsourilo, 2000), 

3 

Regression line 


Range of observed weekly rates from 2007 to 2009 for: 1) two std-MDMs at 40 g/cm 2 , 2) six std 
MDMs at 10 g/cm 2 , 3) two enhanced MDMs at 40 g/cm 2 , 4) one enhance MDM at 10g/cm 2 


* SV aspect ratio = Sensitive volume aspect ratio = T/W 


Table 2: HCOR CMOS SDRAM upset rates are in reasonable 
agreement with FOM predictions 


Device Test Data 

Number 
of bits 

per 

HCOR 

Median 
Shieldin 
g Mass 

Device 

FOM 

FOM 

Calculate 
d Upsets 
per week 

FLUKA Sphere Target 
RCC calculated 
upsets per week per 
MDM 

(SV aspect ratio)* 

Range of Observed 

In-Flight 

Upsets per week 

Samsung 

KM44S32030T-GL 

(high 






















Shuttle 

Flights 

Device 

Median 
shieidin 
g Mass 

Observation 

Time 

Device 

FOM 

Weibull 















15.30 1 6ED2B52 1 
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Table 6 


Vehicle 

System 

Device 

Observation Period 

Median Shielding Mass 

Cassini (interplanetary near 1 
AU) 

Solid state 
recorder 

OKI (4Mx1) DRAM 
640/SSR 

March 2000 

3.4 g/cm 2 

(NASA/JPL estimate) 

SOHO (Interplanetary near 1 AU) 

Solid state 
recorder 

Tl (4Mx1) DRAM 
SMJ44100 

1996 ±2001 

1.0 g/cm 2 

(ESA estimate) 

SOHO (Interplanetary near 1 AU) 

GOLF instrument 

ATMEL CP65 656EV-45 
32kx8 SRAM 

1996-2001 

1.0 g/cm 2 

(ESA estimate) 

Mercury Messenger 
(Interplanetary near 1 AU) 

Solid state 
recorder 

CMOS SRAM 
(4M x 1) 

(ASIC SEE/RAD 
KDUG 3 65$onn'n 

2004 to 2006 

1.0 g/cm 2 

(NASA estimate) 

ETS-V(GEO) 

Technical Data 
Acquisition 
Instrument (TEDA) 

NEC (64k xl) CMOS 
SRAM 

p PD4464D-20 

1987-1997 
2345 days sol max 
1230 days sol min 

5.8 g/cm 2 

(JAXA estimate) 

Thuraya (GEO) 

Digital Signal 
Processor (DSP) 

0.25 U SRAM, 
ASIC, IBM SA-12 
Library 

2001 - 2007 

0.68 g/cm 2 

(Boeing Satellite Development Center 
estimate) 


Vehicle 

Device 

Median Shielding 
Mass 

FOM Caiculated 
Upsets 
per bit day 

FLUKA 

(RCC or isotropic (iSO)) 
Caicuiated 
Upsets per bit day 
(SV aspect ratio) 

In-Flight 
Upsets 
per bit day 

References 

Thuraya (GEO) 

ASiC DSP 0.25 










































54 


GEO/interpianetary spacecraft device weibull 
parameters ' 


Thuraya DSP 0.25 micron CMOS megagate ASIC 

cross section vs LET, aspect ratio, T/W = 0.7ttansen et al 2007 


lodsp: = 2.7 adspsat := 6.3 10 


wdsp : = 20.6 zdsp:=1.2 


odsp (1) : = Gdspsat • 


1 - exp ■ 


1 0.1 - lodsp 
wdsp 



ET5-V. 64K SRAMNECp PD4464D-20; 8 devices in ETS-V 
SEL masked SEU in ground tests - SV depth 1 micron 
aspect ratio = T/W = 0.05, Goka, 1 998 

S: = 0.5 cxsatETS: = -^- wETS:=15 zETS: = 2.9 

64000 


aET S(l) : = asatET S- 1 



l-loETS ^I 
wETS J J] 


* Note ± Weibulf parameters are either taken directly 
from references or were produced by nonlinear least 
squares fitting to heavy ion test data published in the 
references. The Pearson R correlation coefficients are 
greater than 0.9 in all cases. 


Mercury Messenger, aspect ratio approximately 1 , Reed et. al. 2007 


lomm: = 0.3 Gsatmm : = 4* 1 0 


8 


wmm: = 60 zmm: = 6 


omm : = crsatmm- 
i 


1 - exp - 


( LETmm- lomm 
i 


NZimJ] 


V wmm j 


il 


Cassini Solid State Recorder, isotropic target - Swift 200 


Iocs := 0.5 creates :=410 


wcs : = 32 zcs := 3 


ores . : = creates 

i 


1 — exp - 


( LETcs - Iocs 
i 


x ZCS "| 


wcs 


il 


SOHO T1SMX441 00-80 4 M x 1, thinckess = 2 microns; Ha rboe -Sorensen, 2002 


loti: =0.7 ersatti : = 5- 1 0 


wti : = 1 5 zti: = 2.7 


oti. : = cisatti- 


1 - exp - 


LETti - loti ^ 
1 


2ti ] 


wt* y 


iJ 


SOHO MHS CP65656EV 32Kx 8 SRAM; thickness = 2 microns - Ha rboe -Sorensen, 2 


locp:=1.9 asatcp:=610 


7 


acp . := osatep 


1 - exp 


wep : = 17 zep := 1.2 
^zep ] 

IJ 


^LETcp - locp 
L V 


wep J 
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